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BOUNDARY VALUE PROBLEMS ON A HALF SIERPINSKI 

GASKET 

WEILIN LI AND ROBERT S. STRICHARTZ 



Abstract. We study boundary value problems for the Laplacian on a domain 
Q consisting of the left half of the Sierpinski Gasket (SG), whose boundary 
is essentially a countable set of points X. For harmonic functions we give 
an explicit Poisson integral formula to recover the function from its boundary 
values, and characterize those that correspond to functions of finite energy. 
We give an explicit Dirichlet to Neumann map and show that it is invertible. 
We give an explicit description of the Dirichlet to Neumann spectra of the 
Laplacian with an exact count of the dimensions of eigenspaces. We compute 
the exact trace spaces on X of the L 2 and _L°° domains of the Laplacian on 
SG. In terms of the these trace spaces, we characterize the functions in the 
L 2 and L°° domains of the Laplacian on CI that extend to the corresponding 
domains on SG, and give an explicit linear extension operator in terms of 
piecewise biharmonic functions. 
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1. Introduction 

The Laplacian on the Sierpinski Gasket was first constructed as a generator of 
a stochastic process, analogous to Brownian motion, by Kusuoka [6] and Goldstein 
[3j . An analytic method of constructing the Laplacian on the Sierpinski Gasket as 
a renormalizcd limit of graph Laplacians was later developed by Kigami [4] . With a 
well defined Laplacian, it is possible to study differential equations on the Sierpinski 
Gasket, although strictly speaking, these are not differential equations. 

Harmonic functions on the Sierpinski Gasket have been studied in detail and 
the Dirichlet problem on the entire gasket reduces to solving systems of linear 
equations and multiplying matrices. However, there has been little research into 
boundary value problems on bounded subsets of fractals, except for [8], [9] and 
[13] that consider domains generated by horizontal cuts of the gasket. Hence we 
believe it is appropriate to begin our exploration by studying the Dirichlet problem 
on a boundary generated by a vertical cut along one of the symmetry lines of the 
gasket. This is the simplest example of a boundary given as a level set of a harmonic 
function. We hope our results give insight into more general techniques for solving 
the Dirichlet problem and other boundary value problems on more general domains. 

Most of our results are applications of Kigami's harmonic calculus on fractals 
to our half gasket. His theory includes many mathematical objects specific to the 
world of fractal analysis, such as renormalized graph energies, normal derivatives 
and renormalized graph Laplacians. We will present some notation as we proceed, 
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but for precise definitions and known facts (in particular the results that we call 
Proposition), see textbooks [5] and [TTj . 

The Sierpinski Gasket, denoted SG, is the unique nonempty compact set satisfy- 
ing SG — U ? =o FjSG where Fj are contractive mappings given by FjX = (x + qj)/2 
and qj are the vertices of an equilateral triangle. Following convention, the bound- 
ary of SG is defined to be Vo = {loi 1i,Q2}- Hence boundary in our language differs 
from the standard topological definition of boundary. Using the mappings Fj , we 
can iteratively generate a set of vertices V m where m depends on the number of 
times we apply Fj. From V m , we can find a graph approximation r m . See Figure 
11.11 for an illustration. Notice how the boundary points {qj} are oriented and we 
keep this orientation for the entire paper. 
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Figure 1.1. Left to right: r ,ri,r 2 ,r 3 of SG 

We work on the domain 51, which can be defined in terms of the level sets of a 
harmonic function. Let h s be the skew symmetric harmonic function with boundary 
values {h a (q ),h s (q 1 ),h a (q 2 )) = (0,1,-1). Then 51 = {x e SG\V :_h s (x) > 0} 
and dn = q U q x U X where X = {x € SG\ V : h s (x) = 0}. We write fi = Ou90. 

Figure 11.21 provides an illustration of 51, which is precisely the left half of SG 
including the points on the symmetry line. In the figure, we labeled the points 
x m = F™~ l F 2 qi and y m = F™qi. Note that X = {x m }™ =1 , so each x m is 
important for obvious reasons. Each y m is important topologically because the 
removal of any y m turns 51 into a disconnected set. 

We also labeled the open sets Y m = FS % ~ 1 F 1 (SG\ V ). Note that dY m = 
{x-m, J/m-ij Vm} and we write Y m — Y m U dY m . Y rn is classified as a cell because 
a cell is defined to be the image of SG under any compositions of contractive 
mappings Fj. Thus 51 = (J TO=1 Y m , which is an almost disjoint union. 

Although 51 is not globally self-similar because 51 cannot be written as a union 
of smaller copies of itself, it is locally self-similar because each Y m is a fractal. The 
retention of this local property is extremely important for our analysis because any 
result regarding SG also holds for Y m with a proper normalization factor. 
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Figure 1.2. A decomposition of Q 

In the later sections, we will be interested in restriction and extension operators. 
Hence, we need to label points on the other half of the gasket. Let z m and Z m 
the reflections of y m and Y m respectively across the symmetry line containing X. 
Then SG = {J m (Y m U Z m ) is an almost disjoint union and this decomposition will 
be useful in the later sections. 

We begin by studying the Dirichlet problem on $7: 

Am = on ft, 

(1.1) ^ u(qi) = ao on dfl, 

u(x m ) = a m on dfl, 

where A denotes the (Kigami) Laplacian with respect to the standard measure, 
u : tt — >• K is the unknown, and {a m }m=o ^ s tne boundary data. Notice that we 
do not prescribe boundary data at qo even though qo £ dQ. This is by preference 
and is inconsequential because for almost the entire paper, we will assume {a m } 
converges. We will refer to (jl.ll) as the BVP. 

In Section [2] we construct a solution to the BVP using the harmonic extension 
algorithm, which we explain in that section. The space of C(f2) solutions to the 
BVP is one-dimensional, but in general, the solutions blow up at go- We show that 
if the boundary data converges, then we can find a C(fi) solution that is unique in 
this function space. 

In Section [3l we study the graph energy of the C(fi) solution to the BVP. Al- 
though its energy is complicated, the culminating theorem presents an equivalence 
between finite energy and the normalized summability of the the boundary data. 
In fact, finiteness depends only on how quickly the data converges and not on the 
limiting value. 

In Section |4l we show that given stronger assumptions on the boundary data, 
we can obtain the existence of normal derivatives on dQ. In particular, we are 
interested in the behavior of the normal derivatives on X . The normal derivatives 
of the C(fl) solution on X can be found in terms of the boundary data. This 
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relationship allows us to define a Dirichlet to Neumann map and we show that this 
map is invcrtible. 

In Section [5l we discuss both Dirichlet and Neumann eigenfunctions on f2. For 
more information on eigenvalues and eigenfunctions on fractals, see [2] and [10] . 
There are no new eigenfunctions on f2, but for a fixed eigenvalue, its multiplicity 
on Q is different from its multiplicity on SG. For each eigenf unction, we count the 
dimension of its eigenspace. 

Section [6] and Section [7] are closely related to each other. We define a restriction 
operator that maps a function to its restriction to and normal derivatives on X. 
We characterize the function spaces dom L 2 A (SG) and dom^oo A (S'G) in terms of 
the restriction operator. Using this result, we provide necessary and sufficient 
conditions for extending functions in dom i 2A(J7) and dom£=oA(f2) to biharmonic 
functions in dom i 2A(5'G) and dom^oo A(SG) respectively. 

Section [5] acts as an appendix and in this section, we prove numerous lemmas 
about Green's functions and special types of sequences and series. Since these 
results are used in multiple sections and are purely technical lemmas, we have 
decided to place them in its own section. While the sequence and series lemmas 
may not be new, we have not found them in previously published work. 



2. Solution to the Boundary Value Problem 

We begin this section by discussing the graph energy. The energy plays a central 
role in fractal analysis on SG because other objects such as harmonic functions, 
normal derivatives and the Laplacian, are defined in terms of the energy. Given 
a fixed value of m and a real valued function u on SG, the (renormalized) graph 
energy of level m is 

£ m (u)= J2 (lynx)-u(y)} 2 , 

where x ~ m y means x and y are in the same cell of level m. The graph energy of 
u is £(u) = lim m _>. 00 £ m (u), allowing the value +oo. 

Given boundary conditions, we define a harmonic function to be the unique 
function that minimizes the graph energy subject these constraints. Additionally, 
our suggestive use of the word "harmonic" is justified: harmonic functions as min- 
imizers of energy are equivalent to functions that satisfy the differential equation 
Am = 0. The Laplacian A is defined in Section @] 

The simplest tool for constructing harmonic functions subject to boundary con- 
ditions is the harmonic extension algorithm. For a function u defined on V m , we can 
define its harmonic extension to V m +\ as follows. Let {vj} be the three boundary 
points of a cell with {u(vj)} given. Then the harmonic extension of u to the three 
new points is shown in Figure [2T] It is not difficult to see that given u on V m , this 
is the unique extension that minimizes the graph energy at level m + 1 . 

We can apply the harmonic extension algorithm infinitely many times and the 
resulting function on V* = [J m V m will be harmonic. It is not difficult to see 
that functions generated by the harmonic extension algorithm must be continuous. 
Furthermore, V* is dense in SG and so for continuous functions, it suffices to define 
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them on a dense subset. Thus, we say a harmonic function is determined by its 
boundary values. 

u(vq) 




u(vq)-\-2u(v\)+2u(v2) 
5 



Figure 2.1. Harmonic Extension Algorithm 

We can use the harmonic extension algorithm to construct a solution to the 
BVP. Any harmonic function on Y m is determined by its values on dY m . Since 
= |J m Y m , any harmonic function on Q is determined by its value at the points 
{x m } and {y m }- In the following lemma, we see that there are additional constraints 
we must take into account. 

Lemma 2.1. Fix m > 2. Let u be a continuous piecewise harmonic function with 
boundary data given by hl.l\) . Then Au(y m ) = if and only if 

(2.1) u(y m ) = —u{y m -\) - -u{y m - 2 ) - a m - -a m -i- 

o 

Proof. Consider the level m approximation of Y m ^\ U Y m . The value of u at the 
midpoint of y m ~\ and y m -2 and the midpoint of y m ~i and x m -\ are determined by 
the harmonic extension algorithm, shown in Figure [2~121 If Au(y m _i) = 0, then u 
satisfies the mean value property at y m -i- Thus, u{y m -i) is the average of its four 
neighboring points in V m and simplifying that equation yields (|2 . 1 [) . Conversely, if 
(|2.ip holds, then it is straightforward to check that Aw(?/ m _i) = 0. □ 



Theorem 2.2. For every choice of convergent boundary data {a m } ; there is a one 
dimensional space o/C(f2) solutions to the BVP. Given a parameter \, the solution 
to the BVP u\ is the harmonic extension of u\(x m ) — a m , u\{y\) — A and 

(2-2) u x {y m ) = r i F m (X) + ^-G m (X), 



where 



1 / 1 

F m W = j^ I 5A-o -ai - 18j^^afc 



G m {\) = — I -5A+15a + 15ai+4^5 fe a fc 

V fc=2 / 
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2u(y m -i)+2u(y rn - 2 )+a m -i 
5 



u{y m -i) 



u{y m ) 




2u(fc-i)+<i(j/„i-2)+2a m -i 



drn-l 



Figure 2.2. Harmonic extension 



Proof. By Lemma 12.11 u\ must satisfy the recurrence (|2.1[) . The recurrence is 
linear, so we can formulate the equation in terms of matrices. Define 





"I "I. 



and B 



r o 


1] 


3 
. 5 


16 

5 . 



Then the recurrence can be written as 



u\(y m ) 
u\(y m +i) 



B' 



fe=i 



777 — k 



Solving the system, we find that 



u\(y m ) = 3' 



771 — 1 



5A — ao — 



2^ 9k 



Ck 



fe=l 



— (- 

5 m I 14 



flfe+i 



' rn-l \ 

-5A + 15a + ^2 ^ Ck 

v fe=l / 



where Ck = Sa/j+i + 3afc. We want our formula in terms of at rather than Cfc, so 
substituting 

777—1 777 

j^ 5 fc c fc = 4^ 5 k a k + 15oi - 5 m 3a m 

fe=l fc=2 

and 

777 — 1 



f, 1 ^ 1 

E^ c * = 18 E^ 



i 



3 fc -» ~ - ^ 3* °* + ai " ^77 3a ™ 

fe=l fc=2 

into the previous equation for u\(y m ) yields (|2.2[) . Extending these values by the 
harmonic extension algorithm uniquely yields a harmonic function u continuous on 



Since u> is a linear combination of a 3 m term and a l/5 m term, u> may blow up 
at go- Naturally, we ask whether we can find a A such that u\ is continuous on fi. 
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Lemma 2.3. Suppose u\ £ C(O) satisfies the BVP for convergent {a m }- Then 
u\ G C(TT) if and only if 

(2.3) lim u\(y m ) = lim ^(im)- 

m— >oo m— foo 

Proof. Suppose u\ e C(f2) solves the BVP. Then ua is continuous at go, which is 
equivalent to (I2.3[) . Conversely, it is easy to see that go is the only point at which 
U\ can be discontinuous. Then (|2 .3[) implies ua is continuous at go, which shows 
that u x e C(0). □ 

Theorem 2.4. If a m — > as m ^ oo, i/ien i/ie function u given by the harmonic 
extension of ' u(x m ) = a m , 



(2.4) «(s/i) = - L + 01 + 



oo \ 

18 Eg* * 

fc=2 / 



1 / g-1 2 ™ , \ 9^ 1 

^r K - y E 3*°* + -j E 5 afc J + -j E 3T 

V fe=i fe=i / fc=i 



and (for m > 2 j 

(2.5) w(y m ) = ^|ao-^^ + ^ 5 fc a fc I + "- ^J ^a m+fc 

V fe=i fe=i / fc=i 

solves the BVP. Furthermore, this function is the unique solution in C(fi). 

Prw/. Substituting (|2.4|) into (|2.2|) yields (|2.5|l . By triangle inequality, 
1 / 9°°1 2 m \ 9 °° 1 

|«(»m)| < ^r M + 7 E 3* l°*l + y E 5 "l a fcl + 7 E 3* l°™+*l- 

\ fe=i fe=i / fc=i 

We claim that |u(y m )| — > as m — > oo. Clearly the first term tends to zero in the 
limit. The second term tends to zero because convergent sequences are bounded. 
Since both the boundary data and l/5 m converge to zero, for all e > 0, there exists 
M such that for all m > M, we have \a m \ < e and l/5 m < e. For m > M, we see 
that 

oo _. oo _. 

E^i a ™+fei^ £ E^ 



3fe l"™-t-«l — "/__, gfc 2 

fc=l fc=l 



and 

M 



^E 5fc KI = ^rE 5fe KI+ E 5 fc -™K|<C7 l£ ( max : |o fc |)+G a e 



fe=l fc=l fe=Af+l v ~ ~ 

Therefore u satisfies condition (I2.3[) and by Lemma |2~51 u <E C(0). Since harmonic 
functions that are continuous up to the boundary satisfy the maximum principle 
[13] , uniqueness follows from the standard uniqueness argument for linear differen- 
tial equations that satisfy the maximum principle. □ 

Corollary 2.5. If a m — > A as m — > 00 for some constant A, then the function u 
given by the harmonic extension of u(x m ) = a m , 



(2.6) u{yi) = - I a + Oi + 



00 1 \ 

18 EgH 

fe=2 / 



and (for m > 2) 
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1 / 9-1 2 ^ cfc \ 9^ i 

\ k=l fe=l / fe=l 



(2.7) u(y m ) = — I oo - - 2^ ^a k + - ^ 5 a k I + j ^ ^fe a ™+ fc 

\ fc=i fe=i / fe=i 

solves the BVP. Furthermore, this function is the unique solution in C(Tt). 

Proof. Consider the modified BVP 

Am = on f2, 

(2.8) { u( qi ) = oo - A on 9fi, 

u(x m ) = a m — A on <9f2. 



Since a m — A — > 0, the hypotheses of Theorem 12.41 are satisfied. Then there exists 
w G C(O) that solves (|2.8|l and the formula for w(y m ) is given by (|2.5|) under the 
map afc i— > aft — A. By construction, the function u = w + A solves the BVP with 
u G C(£l). The maximum principle implies that u is unique. □ 

3. Energy Estimate 

In this section, we look to answer questions regarding the energy of the C(fi) 
solution to the BVP. In particular, is the energy always finite and if not, can we 
characterize functions of finite energy in terms of a condition on the boundary 
data? Our main theorem shows that harmonic functions on $7 do not necessarily 
have finite energy and provides a simple characterization. 

Given a function u, we say u G dom£ if and only if £ (u) < oo. Following standard 
notation, domo£ is the space of functions that have finite energy and vanish on the 
boundary Vq. It is known that dom£ c C(SG) and in fact, is a dense subset. 

Suppose u is a piecewise harmonic function on f2 that is harmonic on each Y m 
with data given by (|1.1[) . Then the energy of u restricted to Y m is constant after 
level m and is determined by u(y m ), u(y m -i), and a m . It follows that 

/5\ m 
£(")ly m = ( o ) [( u (y™) - u(y m -i)) 2 + (u(y m ) - a m ) 2 + (u(y m -i) - a m ) 2 ] , 

where it is understood that tt(yo) = u (li) = o-o- Then £{u) is the sum of the energy 

of each cell, 

(3.1) 

oo / _ n rn 

£( u ) = X! ( q ) [( u (^™) - u(y m -i)) 2 + (u(y m ) - a m ) 2 + (u(y m -i) - a m ) 2 ] . 

m=l ^ 

If we add the additional assumption that u G C{Vt) solves the BVP, then an equation 
for £(u) as a function of {a m } can be obtained by substituting (12.6[) and (J2.7I) into 
(|3.ip . However, £(u) is series of quadratic terms of series, which is too complicated 
to analyze directly. Instead, we estimate it. 

Lemma 3.1. Suppose u G C(fi) solves the BVP with convergent {a m }. Then we 
have the energy estimate 

Ci y^ I - J (Om+i - a m ) 2 < £(u) < C 2 ^2 ( o ) i a m+i - a m ) 2 < oo. 

m=l ^ m=l ^ 
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Proof. We prove the lower bound first. By ignoring the first term of (|3.1|) . we have 

a m ) 2 + (u(y m -i) - a m ) 2 ] ■ 



/ K\ 



Using basic calculus, we find that u{y m ) = (l/8)(5a m +i + 3a m ) minimizes the 
previous series. Substituting this value of u(y m ) into the previous inequality, we 
obtain 

00 /5\™ 5 5 

m— 1 ^ ^ 

For the upper bound, consider the piecewise harmonic function w given by the 
harmonic extension of w(x m ) = w(y m ) = a m and w(qi) — oq. Since u is a global 
harmonic function while w is a piecewise harmonic function, we have £{u) < £{w). 
Note that £{w) is given by (|3.ip because it; is a piecewise harmonic function satis- 
fying the boundary conditions. Then 

„, n ^ x v 2 ^ AY" 10, N2 10, , 2 

£{u) < £[w) = 2^ I 0) y(am+i-a m ) +y(ai-a J < 00, 

m— 1 ^ ^ 

which completes the proof. □ 

Theorem 3.2. Suppose u G C(J7) solves the BVP with convergent boundary data 
a m — > A. Then u £ dom£ if and only if ||(5/3)" l ' 2 (a m — A)||^2 < 00. Additionally, 
we have the upper bound £(u) < C||(5/3)" l / 2 (a m — ^4)||^. 

Proof. Suppose u e C(fl) solves the BVP with convergent boundary data a m —> A. 
Lemma [3~T1 says that £{u) < 00 if and only if ||(5/3) m/ ' 2 (a m+ i — a m )||#2 < 00. 
Applying Lemma [531 yields the desired statement. □ 

4. Normal Derivatives 

Although the normal derivative and the (standard) Laplacian on SG are defined 
independently, they are closely connected via the Gauss- Green formula. 

For a continuous function u, its normal derivative at qj € Vq, denoted d n u(qj), 
is defined to be 

(4.1) d n u{ qj )= lim (f) [2u(fcO-<i?W)-«(*7Vi)]- 

We say d n u{qj) exists if the above limit exists. In the special case u is harmonic, 
we have the simplified formula 

(4.2) d n u(qj) = 2u(q J ) - wfe-i) - u{q j+ i). 

The formula for the normal derivative of a harmonic function at a boundary point 
of a cell is similar to the above formula, except we require a renormalization factor 
depending on the level. A junction point is a boundary point of two adjacent cells 
of the same level, and the normal derivative with respect to the cells will differ 
by a minus sign. If we need to distinguish between the two normal derivatives at 
a junction point, we use either (<—,—>•), {/*,•/') or (\,\), corresponding to the 
geometrical notion of a normal derivative. 
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Proposition 4.1. Suppose u £ domA. Then at each junction point, the local 
normal derivatives exist and /* d n u + ^/ d n u = 0. This is called the matching 
condition for normal derivatives. 

The Laplacian of a function is defined in terms of its weak formulation. First, we 
define the (symmetric) bilinear form of the energy: given functions u, v and integer 
m, the bilinear form of the energy is 

/5\ m 

s m (u,v)= 2J (3) i u ( x ) - u (y)]l v ( x ) - y (v)]- 

SG has a unique symmetric self-similar probability measure that we denote dx. 
Then the Laplacian can be defined as follows. Suppose u £ dom£ and / is contin- 
uous. Then we say u £ domA with Au = / if 

£(u,v) = — / f(x)v(x) dx 
Jsg 

for all v £ dom £ (functions in dom£ vanishing on Vq). Since £(u, v) = £(v, u), sub- 
tracting the Gauss-Green formula from its transposed version yields the symmetric 
Gauss- Green formula 

(4.3) / (Auv — uAv) dx — y^^vdnU — ud n v) — 0. 

JsG Vo 

The following result relates the normal derivatives of a function with its Laplacian. 

Proposition 4.2 (Gauss-Green). Suppose u £ domA. Then d n u exists on Vq and 
the Gauss- Green formula, 



£(u, v) = — / Auv dx + y^ vd n u, 
Jsg 



holds for all v £ dom£ . 



For the remainder of this section, we assume u £ C(f2) solves the BVP with 
convergent boundary data. Naturally, we are interested in analyzing the behavior 
of d n u(x) for x £ dil. For all points in SI except qo, the formulas for the normal 
derivatives of u are given by (|4.2[) . Using this equation, with the appropriate 
normalization factor, the normal derivative of u at y m with respect to the cell Y m 
is 

/5\ m 

(4.4) /• d n u(y m ) = I - 1 [2u(y m ) - u(j/ m _i) - a m ]. 

Similarly, the normal derivative of u at x m with respect to Y m is 

(4.5) -> d n u(x m ) = ( - J [2a m - u(y m ) - u{y m -x)}- 

However (|4.2[) does not give us the equation for f d n u(qo) because u is only defined 
on Q. But we can define d n u(qo) in a natural way. 

Lemma 4.3. If u £ domA(SG) , then 

(4.6) tM?o)=2- lim / d n u{y m ). 

rn— too 
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Proof. Write u = u s + u a where u s and u a are the symmetric and skew-symmetric 
parts of u respectively. Since u a \ F m, SG ^ = 0(l/5 m ), we have 

t d n u a {q ) = 2 • lim /■ d n u a {y m ) = 0. 

m— >oo 

For the symmetric part, consider the triangle T m with boundary points {qo, y mi z m } 
and the harmonic function v on T m with v(qq) — v(y m ) — v(z m ) — 1. Applying the 
symmetric Gauss-Green formula (|4.3[) for u s and v, we find that 



4- d n u s (q ) + /■ d n u s (y m ) + \ d n u s {z TI 



Au s dx. 



Notice that /• d n u s (y m ) = \ d n u s (z m ) by symmetry. Using the normal derivative 
matching condition of u at qq, we see that f d n u s (y m ) = — \. d n u s (qo). Making 
these substitutions and taking the limit m — > oo, we find that 

2 • lim /• d n u s (y m ) - t d n u s (q ) = lim / Au s dx = 0, 



because Au is bounded and the measure of T m tends to zero in the limit. 



□ 



Motivated by this lemma, we define j" d n u(qo) for u defined on f2 by (|4.6p . In 
the special case that u G C(Jl) solves the BVP with convergent data, then 



(4.7) 



t d n u(q ) = lim 



V 7 fe=m+l v 7 fc=l 



"A: 



which we obtained by substituting (12.71) into the definition of f d n u(q ). 

Notice that (|4.2[) implies that for harmonic functions defined on SG, its normal 
derivative exists everywhere. However, this is not true for harmonic functions on 
n because the limit in (|4.7[) may not exist. It is straightforward to see that if 
a m = Ai + A 2 (3/5) m + o((3/5) m ), then the limit in (|4T7| exists and equals a 
constant times Ai . It is not clear whether or not the converse of this statement 
holds. 

To find the normal derivatives on X in terms of the boundary data, we substitute 
(|2~7) into (1431) . which yields 



(4.8) 



12 ^ 1 \ 1 / 12^ rJ . 54 ^ 1 \ 

fc=i / \ fc=i fc=i / 



where r\ m — — > d n u(x m ). We can think of (|4.8[) as a Dirichlet to Neumann map 
on X because it maps the Dirichlet boundary data to the corresponding normal 
derivatives. Define the infinite vectors 



Vi 



<Zi 



■1/3" 



n = 



Vi 



a = 



and ao = 6ao 



1/3* 
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and the infinite matrices L = Diag[(5/3) 1 ] and K with entries 



K itj 



Then (14.81) can be written as 



4 5* 8 5« 3J " ' -^ ■>■ 



16 

r7 = yL(I-K)a + a . 

Since we assumed {a m } converges and u £ C(fi), we see that {a m }, {u(y m )} £ i°° . 
Then (|4.5|l implies ||(3/5) m 77 m ||^oo < oo. For this reason, for a real number r, we 
define the space 

^ 00 = {{c m }:||r m Cm ||,o <oo}. 
Then we define the Dirichlet to Neumann map D^ : £°° — > £ 3 / 5 '°° given by 

£> A ra=yL(I-K)a + a . 

Theorem 4.4. The Dirichlet to Neumann map is invertible. 

Proof. We see that D N is a composition of L : f 30 -> ^ 3 / 5 >°° with I - K : i°° -)• ^°° 
plus a translation. The translation is not important and obviously L is invertible 
because it is diagonal. 

It is well known that I — K is invertible if and only if p(K) < 1, where p(K) is the 
spectral radius of K. The sum of the entries of the i-th row is 




Consequently, 



IKI 



V^ v 7 3/^1 ^ i \ 



Since K is a positive matrix, the Perron-Frobenius Theorem for positive matrices 
states that p(K) < ||!EC]| OQ . Thus, p(K) < 1, which shows that I— K is invertible. □ 

5. ElGENFUNCTIONS 

The exact spectral asymptotics on the whole gasket and the structure of the 
spectrum has been analyzed previously [T2]. Motivated by that result, we discuss 
eigenvalues and eigenfunctions on the half gasket. Observe that: 

(1) A Dirichlet eigenfunction on il extends by odd reflection to a Dirichlet eigen- 
function on SG and conversely. 

(2) A Neumann eigenfunction on Q extends by even reflection to a Neumann eigen- 
function on SG and conversely. 
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Thus there are no new eigenvalues on fl because odd eigenfunctions on SG are 
Dirichlet eigenfunctions on $7 and even eigenfunctions on SG are Neumann eigen- 
functions on n. Hence we count the number of even and odd eigenfunctions on 
SG. 

On SG, there are #V m = (3 m+1 + 3)/2 vertices on level m, of which to + 1 lie 
on go U X and three are boundary points Vq. The eigenfunctions with eigenvalue 
A < Co5 m for a specific choice of Co are born on level k < m and are in one-to-one 
correspondence with the graph eigenfunctions on V m , so there are (3 m+1 + 3)/2 
Neumann eigenfunctions and (3 m+1 — 3)/2 Dirichlet eigenfunctions. Thus on SI, 

1 /3 m+1 + 3 
^{Neumann eigenfunctions with A < Co5 m } = - I \- to + 1 

1 /3 m +! — 3 
^{Dirichlet eigenfunctions with A < Cq5" 1 } = — I 

because the to + 1 vertices on q U X contribute even functions to the Neumann 
count while the to vertices on X do not contribute to the Dirichlet count. Note 
that the correction terms to + 1 and —to are of the order log 5 m . This is consistent 
with the observation that <9S1 is zero dimensional. We can be more specific about 
individual multiplicities of eigenvalues on 51. For a set U, define the functions 

N(U) = ^{Neumann eigenfunctions on U}, 

D(U) = ^{Dirichlet eigenfunctions on U}. 

(1) 0-series (constant eigenfunctions) have multiplicity N(il) = 1 and D(S1) = 0. 

(2) 2-series only show up in the Dirichlet spectrum on SG, but they are all even 
so they are absent from the Dirichlet spectrum of SI. Thus, iV(fi) = and 
D(O) = 0. 

(3) 3-series are entirely Neumann eigenfunctions on SG that are born on level 
with multiplicity 2. Then JV(fJ) = 1 and D(Q) = 0. 

(4) 5-series are born on level k where k > 1 for Dirichlet eigenfunctions and k > 2 
for Neumann eigenfunctions. If Sk denotes the number of cycles of level less 
than k, then on SG, we find that N(SG) = S k and D{SG) = S k + 2. For a 
cycle that lies on X, the eigenfunction is odd, so that contributes to D(SG) 
but not to N(SG). See Figure liOl for an example of such a function. Note 
that any unlabeled point means the function is defined to be zero at that point. 
Additionally, of the two extra Dirichlet eigenfunctions on SG, exactly one is 
odd, as shown in Figure IQ1 

The number of cycles of level n is 3™ _1 and exactly one of these lies on X. So 
there are (3 n_1 + 1)/2 odd eigenfunctions and (3 n_1 — l)/2 even eigenfunctions. 
Thus 



N(n) = ^2l-(3 n - 1 -i 



71 = 1 



l d k - 1 
2 V ~ 2 



k 



and 



k 



£>(fi) = ^i(3- 1 + l) + l = ^ 

n—l 

(5) 6-series on SG are born on level k where k > 1 for Neumann eigenfunctions 
and k > 2 for Dirichlet eigenfunctions. We know that N(SG) — #Vfc_i and 
D(SG) — #Vfc_i — 3. Neumann eigenfunctions are obtained by giving arbitrary 
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Figure 5.1. Odd cigcnfunction on F 2 




Figure 5.2. Another odd eigenfunction on F 2 



values on the points in Vk-i, while Dirichlet eigenfunctions are obtained by 
giving arbitrary values on the points Vk-i \Vq. 

To find the multiplicities on 0, we just have to count the even eigenfunctions 
and the odd eigenfunctions. Hence 

1 /Qfc j_ Q \ 1 / rtk _ o 

N(Q) = - ^-^ + k) and D(fi) = - ( — ^- - /,- - 1 



6. Trace Theorem 

Consider the restriction map R given by Ru = {(u(x m ), d n u(x m ))}, where u is 
some function defined on some set containing X . That is, R maps u to its function 
values on X and its normal derivatives on X. In this section, we determine the 
image of dom/,2 A (SG) and dom^oo A(SG) under R. We say that u € dom^A^G) 
if u is continuous on SG and Au <G L 2 (SG), and analogously for u e dom^oc A(SG). 

To simplify notation, we define the following spaces. Define the Lipschitz space 

Lip = {{c m }: there exists M such that |c TO +i — c m \ < M for all m}. 

The norm on Lip/Constants is ||c m ||Lip = inf M where the infimum is taken over 
all M satisfying the previous condition. It follows directly from the definition of 
Lip that {c m } € Lip if and only if there exists M such that \c m — c n \ < M\m — n\ 
for all m and n. 
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We define the following trace spaces: 

Too - {{(a m ,f] m )}: a m = A 1 + A 2 (3/5) m + a' m , ||5 m c4IU~, ||3 m 77 m || Lip < oo}, 

T 2 = {{(a m , Vm )}:a m = A 1 +A 2 (3/5) m + a' m , \\(25/3) m ^a' m \\ e 2, \\3 m ^\ n \\ e 2 < oo}, 

with their respective norms 

\\{(a m , Vm )}\\ Tao = |^i| + \A 2 \ + \\5 m a' m \U~ + 113^11^, 

\\{(a m ,V m )}\\ 2 T2 = l^il 2 + \A 2 \ 2 + \\(25/3r/*a' m \\% + \\3 m / 2 Vm \\%. 

Clearly both trace norms satisfy the triangle inequality. Note that the defined norm 
|| • \\j- 2 makes T 2 a Hilbert Space with the obvious inner product. Similarly, we define 
norms on dom^oo A (SG) and dom^A^G) by 

IM|do mi/ ooA(SG) = IM|l~(SG) + I|Au||loo (SG ), 
ll M lldo mi/ 2A(SG) = IMIl 2 (SG) + \\^ u \\l 2 (SG)- 

In the above definition, we could have replaced || • \\ 2 L2 term with || • \\ 2 L001 but that 
would not be a Hilbert Space norm. 

As suggested by the notation, our goal is to prove that R maps dom£oo A(S'G) 
and dom L 2 A (SG) to their corresponding trace spaces. In Section [71 we will show 
that the mapping is onto. 

Theorem 6.1 (Trace Theorem). 

(1) The restriction operator R: domL<x>A(SG) — » Too is bounded and 

||-Ru||Toc < Gi||u|| L =o (SG) + C 2 ||Au|| L oo (SG) . 

(2) The restriction operator R: dom^A^SG) — > 72 is bounded and 

\\Ru\\r 2 < CiIMIl-(sg) + C* 2 ||Au|| L 2 (SG) . 

The proof of the theorem is technical and rather long, so we split the proof into 
multiple lemmas. Our primary tool will be the Green's formula. Given any function 
u on SG for which Au exists, we can write 



(6.1) u(x) = G(x,y)Au(y)dy + h(x), 

JSG 

where G(x,y) is the Green's function (the definition is given in Section l8.ip and 
h is the harmonic function with boundary conditions h\ v = u\ v . We will use 
the Green's function to relate an arbitrary function to its restriction to X and its 
normal derivatives on X. The derivations are digressive, so we have placed these 
computations into their own section. The important formulas and inequalities are 
given by (J8.4I) . (|8.5p . and (|8.7[) . Note that the definition of the function ^ m is given 
in flOj) . 

Since it is easy to check the conditions for the harmonic function h in (|6.ip . let 
us do that first. 

Lemma 6.2. // h is harmonic, then Rh € %>o and Rh £ 7~ 2 with 

(6.2) \\Rh\\ Too = \u(qo)\ + -|«(«i) + u{q 2 ) - 2u(«b)| + \H<li) - «(<fc)|, 

(6.3) ||ita|| T2 = \u(q )\ + ^\u(q 1 )+u(q 2 )-2u(q )\ + ^=\u(q 1 )-u(q 2 )\. 
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Proof. If h is harmonic, then h is a linear combination of the constant function, 
the skew-symmetric harmonic function (with respect to X) and the symmetric 
harmonic function (with respect to X). Then 

r u(q )\ /1\ M / (^ 

u(gi) \=A 1 1 \+A 2 1 \ +A 3 1 

v«(»)/ W W V ly 

where the coefficients are the coefficients A\, A 2 , and A3 are the weights of the 
constant, symmetric and skew-symmetric functions respectively. Solving the system 
for Ax, A2, A3 in terms of u\ v , we find 

M = u(Qo)> A 2 = ~{u{qi) + u(q 2 ) - 2u(q )), and A 3 = -(u(q{) - u(q 2 )). 
On X, we see that 

(1) a constant function is constant with zero normal derivative. 

(2) a skew-symmetric harmonic function is zero with normal derivative A 3 /3 m . 

(3) a symmetric harmonic function has values A 2 (3/5) m with zero normal de- 
rivative. 

Then h(x m ) =Ai + A 2 (3/5) m and d n h(x m ) = A 3 /3 m . □ 

In the following lemma, we prove the bulk of the dom^co A(SG) case. Proving 
the lemma directly from the Green's formula would be difficult, so we employ the 
following indirect method. For the function values of u £ domino A (SG) on the 
vertical boundary, we prove an intermediate statement about the linear combination 
5u(x m +i)— 3u(x m )- We consider the linear combination 5u(x m +i) — 3u(x m ) because 
the troublesome X)feLi ^fc(l, 2, 2) term of (|8.5|) cancels out in the linear combination 
5G(x m +i, y)~3G(x mi y). Then the intermediate result, coupled with a lemma from 
Section 15721 will give us the desired statement, except for a few estimates which we 
prove without much trouble. 

Likewise, for the normal derivatives of u £ dom^oc A (SG) on the vertical bound- 
ary, we prove an intermediate statement about the linear combination 3n m+ i — r) m 
because the troublesome Ylk=i 3 fe *fc(0, —1,1) term in (|8.7[l disappears in the lin- 
ear combination. The intermediary result, combined with the proper lemma from 
Section [5721 and more bounding, yields the desired normal derivative estimate. 

Lemma 6.3. If u £ dom-L°°A(SG) with u = on Vq, then Ru £ Too and 
(6-4) \\Ru\\ roo <C\\Au\\ L ^ iSG) . 

Proof. Suppose u £ dom^oo A (SG) with Ru = {(a m ,f]rn)}- Using the Green's 
formula (Proposition 18. l| ) on 5a„ l+ i — 3a m and the equation for G(x m ,y) given by 
75)) . after some simplification, we obtain 



5o m+1 - 3a m = i (J\ J [3*^+1(3, l,l)-5* m (-l,l,l)]A«di/ 
Then applying inequality (|8.4[) yields 

|5o m+ i-3o m |<||A«|| £ o.^r|j J |3* m+ i(3,l,l)-5* m (-l,l,l)|d» 



C 

— II 11-^ r m 
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Rearranging the above inequality yields 

||5 m (5a m+ i - 3a m )\\l°° < C||Au||i«. 
Lemma 151)1 implies that a m — A(3/5) m + a' m , where A — lim TO _ i . 00 (5/3) m a m and 

||5"XJ|^o <||5 m (5a m+1 -3a m )||,oo. 
The previous two inequalities immediately yield 

(6.5) \\5 m a' m \\ iao <C\\Au\\ L ^. 

It follows from the Green's formula and standard bounding methods that 

3 J M<U) J sG \G(x m ,y)\\Au\dy<C 1 \\Au\\ L ~ + C 2 \\Au\\ L ~ — . 

Since A = lim m _ > . 00 (5/3) m a m , the above implies that 

(6.6) \A\<C\\Au\\ L °°. 

We use a similar technique to prove the desired statement about the normal deriva- 
tives. Using the equation for r\ m given by (|8.7|) to compute 3?7 TO +i — r) m , we obtain 

3r) m+1 - r\ m = — / [-3* m+ i(5, 1, -1) + 5* m (l, -1, l)]Au dy - 3<p m +i + fm, 
1U Jsg 

where ip m was defined in the lemma. Then 



|3r/ m+1 - r] m \ < C\\ Au|| L oo / |3* m +i(5, 1, -1) - 5* m (l, -1, 1)| dy + \3cp m+1 - ip r . 

JSG 

<C\\Au\\ L °° — , 

— II II -f om ' 

where we used (|8.4j) and (|8.1[) to bound the first and second terms respectively. 
Rearranging, we find that 

p m {3T}m+l-Vm)\U°° <C\\Au\\ L ~. 

The above estimate allows us to apply Lemma 18.71 which gives us 

||3 m 77 m ||Li P H|3 m (37 ?m+1 -77 m )|k=o. 
The previous two inequalities imply 
(6.7) ||3 m ry m || L ip < C\\Au\\ L ~. 

Finally, combining our inequalities (j6.5[) , (J6.6P and (|6.7p . we see that 
\\Ru\\ ToB = \A\ + ||5"Xnll^ + l|3 m r?m||Li P < C||Au||£«,. 
Since a m — A(3/5) m + a' m and ||i?u|| < oo, we conclude that Ru £ Too- O 

In the following lemma, we prove the majority of the dom^A^G) statement 
of the Trace Theorem. We use an indirect approach similar to that of the proof for 
the dom£oo A(S'G) case, except the statements are considerably harder to prove. 
Proving the lemma directly from the Green's formula without proving the inter- 
mediary result would be extremely difficult, mainly because the Cauchy-Schwarz 
inequality is too wasteful for the type of estimate we desire. 

The outline of the proof is similar to that of Lemma [6731 For u € dom^2 A(SG), 
we prove intermediary results about the linear combinations 5a m +2 — 8a m +i + 3a m 
and irj m+ i — 16ry m +i + 5r) m , where as usual, a m = u(x m ) and r\ m = d n u(x m ). These 
linear combinations are written as linear combinations of integrals, but the primary 



18 WEILIN LI AND ROBERT S. STRICHARTZ 

integrand of each linear combination is supported on a set not containing go • This 
support allows us give a more precise estimate, thereby limiting the wastefulness of 
Cauchy-Schwartz. Then applying results from Section [53] and some more bounding 
will give us the desired statements. 

Lemma 6.4. If u £ dom^A^SG) with u = on Vq, then Ru <E T 2 and 

(6.8) \\Ru\\ T2 <C\\Au\\ LHSG) . 

Proof. Suppose u £ dom L 2 A(SG) with Ru — {(a m ,^ m )}. Using the Green's for- 
mula ( Proposition 18. 1[) on 5a m+2 — 8a m+ i + 3a m and the equation for G{x m ,y) 
given by (|8.5j) . after much computation, we obtain 

5a m+2 - 8a m+ i + 3a m = f - J / G m Au dy, 

where we defined 

G m (y) = ^[9* m+2 (3, 1, 1) - 20* m+1 (l, 0,0) + 25tf m (l, -1, -1)]. 
50 

We show that G m is supported on D m = Y m U Y m+ i U Y m+2 U Z m U Z m+ i U Z m+2 . 
Since G m is a linear combination of harmonic splines, we see that G m vanishes on 
Y m i U Z m i for m! < m. Using the harmonic extension algorithm, notice that 

25tf m (l, -1, l){y m+2 ) = 25* m (l, -1, l){z m+2 ) = -9, 
20* m+1 (l,0,0)(y m+2 ) = 20* m+1 (l,0,0)(z m+2 ) = 0, 
9* m+2 (3, 1, l)(z/ m+2 ) - 9* m+2 (3, 1, l){z m+2 ) = 9. 

Thus G m (]Jm+2) = G m {z m +2) = and consequently, G m vanishes on Y m > U Z m > for 
m' > m + 2, which proves that G m is supported on D m . Taking advantage of the 
support of G m , we can write 

5a m+2 - 8a m+ i + 3a m = ( 7 ) / G m Au dy, 



Applying Cauchy-Schwarz and inequality (18.41) on the above equation yields 

|5a„ l+2 - 8a m+1 + 3a m | 2 < C\\ Au\\ 2 L2{Dm) f — 
By definition of D m and linearity of the integral, we have 

m+2 

\\ Au \\h(D m ) = Yl W Au \\ 2 LHY k uz k )i 
(6.9) k=m 



|Am||| 2(SG) = ^||A u 



2 

L2(y fe uz fc )- 
fe=i 



Using the upper bound on |5a m+ 2 — 8a m+ i + 3a m | 2 and the above two equations, 
we obtain 

|!(25/3) m/2 (5a m+ 2 -8a m+ i +3a m )|| P < C|| Au|| L 2 (SG) . 

This estimate allows us to apply Lemma [8.101 Thus a m = A\ + A 2 (3/5) m + a' m , 
where A\ = lim™-^ a m , A 2 — lim m ^ 00 (5/3) m (a m — Ai), and 

||(25/3) m / 2 a^|| £2 < C||(25/3) m / 2 (5a m+2 - 8a m+1 + 3a m )||^. 
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The above two inequalities immediately yield 

(6-10) ||(25/3r/XJL 2 < C\\Au\\ LKS G)- 

We claim that Ai = and \A 2 \ < C||Au|| i 2(5 G ). Applying Cauchy-Schwarz to the 
Green's formula for a m , we find that that 

'5\ m 1 

- j \a m \ < Ci\\Au\\ L 2 (S G) +C 2 ||Au|| L 2 (SG) ^ ;r ^. 

The above inequality implies that Ai = and 

(6.11) \A 2 \<C\\Au\\ L 2 {SG) . 

We use a similar argument to prove the estimate on the normal derivatives. Using 



Lemma 18751 to compute 3r] m +2 — 16?7 m +i + 5r] m , we see that 

3r? m + 2 - 16r/ m+ i + 5ry m = / <f> m Au dy - (3</? m+2 - 16ip m +i + 5<p m ), 

JSG 

where we defined 

$m = -^[-3* m+2 (5, 1, -1) + 10* m+1 (8, 1, -1) - 25* m (l, -1, 1)]. 

We show that $ m has support on D m as well. Since $ m is a linear combination 
of harmonic splines, $ m vanishes on Y m i U Z m i for m! < m. Using the harmonic 
extension algorithm, we have 

- 25* m (l, -1, l)(y m+2 ) = 25* m (l, -1, l){z m+2 ) = 1, 

- 10* m+ i(8, 1, -l)(y m+2 ) = 10* m +i(8, 1, -l)(z m+2 ) = -2, 

- 3* m+2 (5, 1, -l)(y m+1 ) = 3* m+2 (5, 1, -l)(z m+ i) = -3. 

Thus, $ m (j/ m + 2 ) = Qm{z m +2) = and consequently, <& m vanishes on Y m > U Z m > for 
m! > m + 2. Using the compact support of $,„, we can write 

3r? m+2 - lSrjm+i + 5r] m = / $ m Aii dy - (3<^ m+2 - I6ip m+1 + 5<^ m ), 

It is straightforward to find an upper bound on the linear combination of tp m terms. 
Using Cauchy-Schwarz and inequality (|8.2|) . we obtain 

|37j m+2 - 16?7 m+ i + 5?7 m | 2 < C (\<p m+ 2\ 2 + \<f m +i\ 2 + \f m \ 2 ) < C\\ Au\\ 2 L2{Dm) — . 
Using Cauchy-Schwarz and inequality (|8.4p . we find that 

<& m Au dy 



< || M\h( Dm ) J \®mf dy < C\\Au\\l 2(Dm) ^ 



Combining the above two inequalities and (16.91) yields 

(6.12) ||3 m / 2 (3 ??m+2 - 16?7 m+1 + 5^)11,2 < C\\Au\\ L 2 (SG) . 

The hypothesis of Lemma T8. Ill is satisfied, so we have r\ m = 5 m A + Tj' m with 

(6.13) \\3 m/2 v' m \\ < Ci(ri2 ~ 5?7i) 2 + C 2 ||3 m (377 m+2 - 16r /m+1 + 5ry m )||^. 
However, applying Cauchy-Schwarz to (|8.7[) yields 



I 7 ?™! < C ll Au ll ^(SG)^^- 
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This forces A = and so r\ m = n' m . Note that the above bound provides the upper 
bound (772 — 5?7i) 2 < C\\ AuW^tgQ)- Combining this inequality with (|6. 12|) and 
(|iTT5|l yields 

(6-14) ||3™/ 2 ry;j| 2 <C||A U ||2 2(SG) . 

Finally, using (jOO]) . ([6lT|) and (|6T4]) , we see that 

||flu|£ a = l^l 2 + |A 2 | 2 + ||(25/3)™/ 2 a^|| 2 2 + ||3 m /V|| 2 2 < C||A U || 2 2(SG) . 

Since a m = A2(3/5) m + a' m and ||i?u||^- < 00, we conclude that Ru G Ti- □ 

Finally, we have the necessary results to prove the Trace Theorem. 

Proof of the Trace Theorem. Suppose u G dom^oo A(SG) or u G dom^2A(5G), 
and Ru = {(a m , r] m )}). Let h be the harmonic function determined by the boundary 
values h\ v — u\ v . Let w — u— h, and note that Aw = Aw and w = 0on Vq. The 
Green's formula states that 



u(x) = h(x) + / G(x,y)Aw(y) dy. 
Jsg 

(1) Suppose u G doixiL=o A (SG). Using triangle inequality on u = w + h, the 
estimate (|6.2[) applied to h, and the estimate (|6.4j) applied to w, we find that 



pullToo < k('Zo)l + 2l u (9i) + M (92)-2u(go)l+2l u ('?i)- M (92)|+C||A M || LTC(SG) . 

(2) Suppose u G domi2A(5G). Using triangle inequality on u = w + h, (16. 3[) 
applied to /i, and (16.8[) applied to w, we find that 



|^w||r 2 < l - "((7o)| + 2 -|u(gi)+u(g 2 )-2u(qo)| + --^=|u(gi)-'"((72)|+C , ||Au|| L 2 (SG) . 

□ 



7. Extension Operators 

In this section, we present two different extension theorems. The first extension 
will be a right inverse to the restriction map R. The second extension will map 
solutions to differential equations on the half-gasket to a well-behaved function 
on the whole gasket. The ideas behind the two extensions are similar, but with 
different computations and formulas. In order to construct the desired extensions, 
we will require the following result. If will give us the exact conditions under which 
a piecewise function is in the domain of the Laplacian. 

Proposition 7.1 (Gluing Theorem). Let u and f be defined by gluing pieces {uj} 
and {fj} (j = 0, 1, 2), with Auj = fj on FjSG. Then u G domA with Au = f if 
and only if fj(FiQj) = fj(FjQi) (i 7^ j) holds for {uj} and {fj} (so u and f are 
continuous) and the matching conditions on normal derivatives hold at the three 
points. 
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7.1. The Inverse Operator to R. We seek a linear extension operator E that 
is a right inverse of the restriction operator R. The desired extension will satisfy 
E: Too — > domL= A(SG) and £: T2 ^ dom L 2 A(SG). In order to construct this 
extension operator, we study piecewise biharmonic functions. Biharmonic functions 
satisfy the differential equation A 2 u = and in particular, biharmonic functions 
satisfying Au = C for some constant C is a four-dimensional space on SG. One 
way to specify a constant Laplacian function on SG is to specify the value of the 
function on Vq and the constant. 

Lemma 7.2. Suppose An = C on some cell of level m with boundary points 
P0iPiiP2- Then the outward normal derivative of u at pj is 

(7.1) d n u(pj) = f - J [2u(pj) - u(p j+ i) - u(pj-i)] + ^+i- 

Proof. Let v be the harmonic function on the cell with the boundary values v(j)j) = 
1 and v(pj + i) = v(pj—i) = 0. Since v is harmonic on a cell of level to, using (|4.1[) 
with the proper normalization, we have d n v(pj) = 2(5/3) m while d n v(pj+i) = 
d n v(j)j-i) = — (5/3) m . Applying the symmetric Gauss-Green formula (j4.3|l . we 
obtain the desired formula. D 

Lemma 7.3. Given any sequences {a m } and {rj m }, there exist a piecewise bihar- 
monic function u on SG and sequences {C' m } and {C m } such that Ru = {(a m , r) m )}, 
Am = C' m on Y m , Am = C m on Z m , and the normal derivative matching conditions 
hold at {x m }, {y m }, and {z m }. 

Proof. We construct two functions u\ and U2 such that ui{x m ) — a m but d n ui(x m ) — 
0, while U2(x m ) = but d n u(x m ) = r\ m . Then the sum u = U\ + u-i will satisfy 
Ru = {(a m , rjm)}- Of course, we must do this carefully so that u satisfies the other 
claimed properties. 

Consider the symmetric piecewise biharmonic function u\ satisfying ui{x m ) = 
a m , ui(y m ) = ui(z m ) = (l/8)(5a m+ i + 3a m ), and Aui = D' m on Y m U Z m with 

D' m = 5 m ( - J (5a m+ i - 8a m + 3o m -i). 

This information determines u\ on Y m UZ m because as mentioned earlier, a constant 
Laplacian function is determined by its boundary values and the value of its Lapla- 
cian. Consequently, u\ is determined everywhere because SG = {J m (Y m U Z m ). 
Using (|7.ip to compute the normal derivatives of u\ at x m , y m and z m , it is straight- 
forward to check that d n ui(x m ) = and the normal derivative matching conditions 
hold. 

Consider the skew-symmetric piecewise biharmonic function U2 satisfying the 
conditions u 2 {x m ) = 0, u 2 (y m ) = -(l/8)(3/5) m (r? m+ i +rj m ), u 2 (z m ) = -u 2 (y m ), 
Am2 = ~E m on Y m and Au2 = E m on Z m , where 

E m = 3 m f - J (3r? m+ i - 16r? m + h-q m -i)- 

Again, these constraints determine u 2 everywhere on SG. Writing down the normal 
derivatives of u 2 at x m , y m and z m using (|7.ip . we see that d n u 2 {x rn ) = r\ m and 
the normal derivative matching conditions hold. 
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Then the function u = u\ + U2 satisfies u(x m ) — a m , d n u(x m ) = f] m , 

1 1 /3\ m 

u(y m ) = o (5a m+ i + 3a m ) - - f - J (jfo+i + J? m ), 

(7.2) ' ^ ' 

1 1 /3\ m 

u(z m ) = -(5a m+ i + 3a m ) + - ( - J (??m+i + ??m), 

A?i = C^j on Y" m and Au = C m on Z m where 
(7.3) 

c 'm = 5 m (I ) (5a m+ i - 8a m + 3a m _i) - 3 m f M (3r? m+1 - 16^ + 5r? m _i) 



C m = 5™ I - j (5a m+ i - 8a m + 3a m _i) + 3'™ ( - j (3r] m +i - 16?7 m + 5?7 m _i). 

Because normal derivatives add linearly, u satisfies the normal derivative matching 
conditions at x m , y m and z m . D 

As a result of the above lemma, we can define the extension operator E which 
maps two sequences {(a m , rj m )} to the function u given in the lemma. This operator 
is well defined because the process described by the lemma generates exactly one 
function for each pair of sequences. Additionally, it is not difficult to see that E is 
a linear operator. 

Theorem 7.4. There exist a bounded linear extension map E: Tx ~^ domi J x,A(SG) 
and E: T — > dom L 2A(SG) with Ro E = Id. 

Proof. Suppose {(a m ,?7 m )} S %>o and let u = E{(a m ,r) m )}. In order to apply the 
Gluing Theorem, we need to check that u is continuous. It suffices to check for 
continuity at qo because u is clearly continuous everywhere else. In order to show 
that u is continuous at qo, we need to show that linim^oo u(x m ) — linim^oo u(y m ) = 
linim^oo u{z m ). Since {(a m , r/ m )} G Too, we have a m — A 1 + A 2 (3/5) m + a' m with 
l|5"XJU < 00 and ||3 m ?7 m || Lip < 00. Then flL2l) reads 

3 /3\ m 1 ,1 /3\ m 
u(y m ) = A 1 + - f - J A 2 + -(5a m+1 + 3a' m )-- ( - J (r} m +i+Vm), 

3 /3\ m 1 1 /3\ m 

u(z m ) = A 1 + -i-j A 2 + -(5a m+1 + 3a' m ) + - ( -J (v m +i+Vm)- 



Taking the limit m —¥ 00 in the above equations, we see that A\ — X\Yn rn ^,oo u(y m ) = 
linim^oo u(z m ) = lininj^oo a m , which verifies the continuity of u at q Q . Recall that 
Lemma 17.31 tells us that u satisfies the normal derivative matching conditions at 
{x m }, {ym} and {z m }. Thus the hypotheses of the Gluing Theorem are satisfied, so 
the theorem implies that Au is well defined. We need to show that Aw e L°°(SG). 
Observe that (17.31) reads 



f','„ : .V ( ( 5a ^ +1 - 8a^ + 3^0 ~ 3 m (~ ) (%„-, - '<>'/- - -»h, -i). 
C TO = 5 m ( - J (5a^ l+1 - 8aJ„ + 3a4 l _ 1 ) + 3 m f - ) (3r/ m+ i - 16ry m + 5r? m _i). 
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Using Lemma [8771 to obtain an upper bound on the normal derivative terms in C m 
and C' m , we find that 

||Au|| LOO < ||C ro || i<)0 + ||d|^ < M x \\S m a' m \\ t «, + M 2 113-^11^ . 

Therefore, E:Toc^ dom L ~A(SG). 

Suppose {(a m ,r] m )} £ 72 and let u — E{(a m , r] m )}. Again, we need to check that 
u is continuous at go in order to apply the Gluing theorem. By definition of 72 , we 
havea m = A x +A 2 (3/5) m + a' m with ||(25/3) m / 2 a^||f2 < oo and \\3 m / 2 r ]m \\ i 2 < oo. 
Then \a' m \ — > and |r7 m | — > 0. By the same argument for the Too case, u is 
continuous at go, hence continuous everywhere. By Lemma 17.31 u satisfies the 
normal matching conditions at {x m }, {y m } and {z m }. Then Au is well defined by 
the Gluing Theorem. Finally, Au £ L 2 (SG) because 



3C 



IIAullJ. = E — 3 t' — < Mi £ (f ) laj 2 + M 2 £ S-I^l 2 . 

771= 1 m=l ^ ' m— 1 

Therefore, E: T 2 -> dom L2 A (SG). D 

7.2. Extensions of Solutions to Differential Equations on 57. The material 
presented in this section is motivated by the classical theory of extending functions 
with Au £ L p on a nice domain in Euclidean space R™ to functions with the same 
property on R™. We ask: 

(1) Given u £ dom^co A(0), does there exist an extension m £ dom^oc A(S'G)? 

(2) Given u £ dom^2 A (57), does there exist an extension u £ dom L 2 A(5G)? 

We present two motivating examples before we proceed to the main extension re- 
sults. 

Theorem 7.5. If u is a nonconstant harmonic function on $7, then its even reflec- 
tion is not in domA. 

Proof. Suppose, for the purpose of contradiction, even reflection extends to a func- 
tion u £ domA(SG). Then Au = on SG, and 

-► d n u(x m ) = 2u(x m ) - u(y m ) ~ u(y m -i), 

<- d n u(x m ) = 2u(x m ) - u(z m ) - u(z m _i). 

Since u is assumed non-constant, both normal derivatives are nonzero. However, 
— > d n u(x m ) = <— d n u(x m ), so the normal derivative matching condition at x m does 
not hold. Therefore, u <£ domA (SG). D 

Theorem 7.6. Suppose u £ C(fi) solves the BVP with ao = C\ and a m — 

(2/3)(3/5)™ l (Gi + C2) for some constants C\,Ci- Then there exists a harmonic 
extension of u. 

Proof. Consider the harmonic function u on SG determined by the boundary values 
u(qo) = 0, u(qi) — C\ and u(q2) = Ci- Simple computation shows that u(x m ) = 
(2/3)(3/5)" l (Ci + C2). Thus, u = u on 57 and Am = 0, which shows that u is indeed 
a harmonic extension. D 

In special cases, such as the one presented in the previous result, there exists a 
harmonic extension. In general, the desired extension will not be harmonic because 
the space of harmonic functions on SG is a three dimensional space so finding a 
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harmonic extension u of u satisfying the infinite number of conditions u(x m ) = a m 
is unlikely. For that reason, we look for a piecewise biharmonic extension. In fact, 
this motivates our study of piecewise biharmonic functions to begin with. To prove 
the existence of an extension, we need the analogue of Lemma 17.31 

Lemma 7.7. Suppose u £ dom^ A(Q) or u £ dom L 2 A(Q). Then there exist a 
sequence {C m } and a piecewise biharmonic function u on SG satisfying u = u on 
fl, Am = C m on Z m , and the normal derivative matching conditions hold at {x m } 
and {z m }. 

Proof. For convenience, we write a m = u{x m ) and rj m = d n u(x m ). Consider the 
function u = u on Q, 

_ 1 1 /3\ m 

(7.4) u(z m ) = -(5a m+ i + 3a m ) + - f - J (r} m +i + Vm), 

and As = C m on Z m where 

(7.5) C m = 5 m ( - J (5a m+ i - 8a m + 3a m _i) + 3 m ( - J (3r) m+1 - 16?y m + 5?? m _i). 

For the same reason as before, these constraints completely determine u on Z m - 
Hence we have defined a function u on SG. 

We claim that the normal matching conditions hold at x m and z m . Using (|7.ip , 

/ 5 \ m C 

•s- d n u{x m ) = I - I [2u(x m ) - u(z m ) - u{Zm-l)} + 



' 5 \ m C„, 

d n u(Zm) = ( - ) [2u(z m ) - u(z m -l) - u{x m )] + 



^ , L ^ v ~ m , ^y^m-ij ^y^mjj i g TO+1 > 

/5\ m+1 C 

\ d n u(z m ) = ( - J [2w(z m ) - u(z m+ i) - u(x m+ i)] + ™ +a - 

It is straightforward to check that our formulas for u(x m ), u(z m ), and C m imply 
the matching conditions hold at {x m } and {z m }. D 

The lemma allows us to define an extension operator. Let Eq be the extension 
operator that maps a function u £ domino A ($7) or u £ dom^Aifl) to the function 
E^u on SG as given in the lemma. This operator is well defined because for each 
u, there is exactly one Enu. It is clear that Eq is linear and that Equ is continuous 
except possibly at qo- 

Theorem 7.8. Suppose u £ domj,^ A(J7) . If Ru £ Too, then Equ £ domL°°A(SG) 
and 

\\A(Enu)\\ L oo (SG) < ||Au|U.o(n) + C r ||.Ru||r 0O . 
The Trace Theorem implies the converse: if Equ £ dom^ A(SG), then Ru £ Too- 

Proof. Suppose u £ dom£ooA(£7) and Ru = {(a m ,^ m )} £ T^. By definition of 77», 
we have a m = A 1 +A 2 {3/5) m + a' m with ||5 m a',J^ < oo and ||3 m ?7 m ||Lip < 00. We 
need to check that Eq is continuous at qo. Observe that (|7.4I) becomes 

/3\ m 1 1 /3\ m 

Enu(zm) = A t + A 2 f - J +g(5a' m+1 + 3a' m ) + - ( -J (vm+i + Vm)- 
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Taking the limit in the above equation, we see that A\ — lim m ^.oo a m — lim m ^oo Enu(z m ). 
This proves that Equ is continuous. By Lemma 17.71 the matching conditions for u 
at {x m } and {z m } are satisfied. This allows us to apply the Gluing Theorem, and 
so A(Equ) exists. 

To prove that Equ G dom£ooA( 1 S , G), observe that 

||5 m (5a„ l+1 - 8a m + 3a m ^)\\ eao < 16 ||5 m a^||^ 
and by Lemma T8.7[ 

||3 m (3r? m+1 - 16r) m + 5r? m _i)||^ < 16 \\T\ n \\ Up . 
Using the above inequalities and the equation for C m given by (|7.5p . we find that 

||A(£bu)|U«. ( n/) = max|C m | < M^a'JU- + M 2 \\3 m r) m \\ Up - 
Then by triangle inequality, 

\\A(E a u)\\ L o° (S G) < ||A«|U«. (n) + M 1 ||5 m aUI*~ + M 2 ||3 m 77 m || Lip , 
which completes the proof. □ 

Theorem 7.9. Suppose u G dom L 2 A (ft) . If Ru G Ti, then Equ G dom L 2A(SG) 
and 

\\&(Enu)\\ 2 LHSG) < \\Au\\ 2 L2m +C\\Ru\\ 2 T2 . 

The Trace Theorem implies the converse: if Equ G dom^i A(SG), then Ru G Ti- 

Proof. Suppose u G dom^2A(rj) and Ru = {(a m ,r) m )} G Ti- By definition of 
75, we know that a m = A x + A 2 (3/5) m + a' m with ||(25/3) m / 2 a;j| £ 2 < oo and 
||3" l / 2 ?7 m ||f2 < co. Then \a' m \ — ¥ and |r? m | —> 0. Using these limits, the same 
argument given in the proof of Theorem 17.81 shows that Equ is continuous. Again, 
Lemma 17.71 guarantees the matching conditions for u at {x m } and {z m } hold. The 
Gluing Theorem implies A(Equ) is well defined. 
To see why Equ G dom L 2 A (SG), we first see that 



<S*J -. WW / O K \ 

\A(E n u)\\l 2(iV) = J2 ^\Cm\ 2 < ^i E T I "' 2 + Ah E 3 ™^ 



771 — 1 m— 1 ^ 771 — 1 

Since \\A(E^u)\\ 2 L2 ^ SG) = ||Ati||| 2 / Q N +\\A(EQu)\\% 2 , n ,y using the above inequality 
gives us 

oo 



II (A^llli.,,,., II. ."Hi,..,,+ U, ^ (|) | a ^| 2 + Af 2 E3" l |^| 2 . 

' 777 — 1 

D 



We can interpret Theorem 17.81 and Theorem 17.91 by the following: Ru G %o 
is the minimal condition for extending an arbitrary function in dom£ooA($7) to a 
function in domiooA(5G) and Ru G 75 is the minimal condition for extending an 
arbitrary function in dom^A^) to a function in dom i 2A(5G). 

A function belonging to dom i 2A(J7) or dom£ooA(£7) is naturally a solution to 
the differential equation Au = f for / G L 2 or / G L°° respectively. Solutions to 
this differential equation can be found using Theorem 18.21 

As a special case of Eq , we can extend harmonic functions u on ft provided that 
Ru G 75 or Ru G Too- Recall that the solution to this differential equation was 
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explicitly given in Section [2l The formula for the extended function will be given 
by (17. 4[) and (|7.5[) . which can be simplified by using the normal derivative formula 
for harmonic functions (|4.2[) and the recurrence relation (j2.ll) . 

8. Appendix 

8.1. Green's Function Formulas. For a given m and a point x E V m \ Vy, let 
ijj™{y) denote the piecewise harmonic spline of level m satisfying ip™(y) = 5 x (y) for 
y E V m and extended harmonically for levels m! > m. Notice that ip™ E domo£ 
because x £ Vo- 

Proposition 8.1 (Green's Formula). On SG, the Dirichlet problem —Am = / on 
SG\Vq and u = on Vo has a unique solution in dom/S for any continuous f, given 
by u(x) = / 5G G(x, y)f(y) dy for the Green's function G(x, y) = lim M ->oo G M (x, y) 
(uniform limit) where 

M 

G M (x, y) = J2 E 9(s, s')i> k 8 (x) ifi(y) 

k=is,s'ev k \v k . 1 

and 

I -(*■* 
ff(s.s') , 

I jq (|) for s, s' E F W K, \w\ = k — 1 and s ^ s' . 

From the Green's formula, we have the following simple observation. 

Theorem 8.2. Let G(x,y) denote the Green's function on SG. Let Ga(x,y) — 
G(x,y) — G(x,Ry) for x,y E f2 where R denotes the reflection. Then Go is the 
Green's function for Q, namely 



JL(§)~ fors = s'£V k \V k _ 1 , 



u(x) = / Ga{x,y)f{y) dy 
Jo. 

solves — Aw = f on Q subject to u\q = 0. 

To simplify notation, we drop the superscript m on functions of the form ip™ , 
ty™ , and ip^ because unless otherwise notated, the superscript index matches the 
subscript index. It follows immediately from the definition that 

2 



(8-1) / IWJ dy = / |V„J dy = / |^J dy = - , . 

JSG JSG JSG 3 

Additionally, since IV'a;™! 2 < IV'xml) we nave 

(8-2) / \iP Xm \ 2 dy=j \ip y j 2 dy=[ \i> z J 2 dy < -JL. 

J SG J SG J SG 

To further simply notation, define the function 

(8-3) <& m (a, 6, c) (y) = a^ Xm (y) + b^ Vm (y) + cip Zm (y) . 

Using (|8.1|) and (|8.2p . we have the estimates 

(8.4) f \^ m (a,b,c)\dy<^ and f |* m (a,fe,c)| 2 dy < §| 

JSG 6 JSG 6 

for constants C\ and C2 depending only on a, b, c. 
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Lemma 8.3. The Green's function evaluated at x m is 

(8.5) G(x m ,y) = - - £tt fc (l,2,2)(i,) + - |) * m (l,-l,-l)(y). 



15 V 5/ ^ v '"" 6 V 5 

v 7 fc=i v 

Proof. Note the following observations: 

(1) If fc > to, then ^J(a; m ) = 0. 

(2) If k = to, then ^^(^m) = 1- If fc = to and s 7^ x m , then i/>™(a; m ) = 0. 

(3) If fc < to with s^t/j. and s 7^ z^,, then ^(s m ) = 0. 

Using these facts, we have 

m— 1 

G(a; m ,y)=^ ^ [fir^, s')i/) yk (x m ) + g{z k , s')ip Zk (x m )\ ip^,(y) 
fc=i S 'evi\v t -i 

s'GU m \u m -i 
Using the harmonic extension algorithm, for k < to, we have 

Vv {Xm) = g f g J and ^ fc ( Xm ) = 3 ( 5 

Since <7(s, s') = if s and s' are in different cells of level fc — 1, we deduce that 

E \9{yk,s')+g{z k ,s')}^ k Ay) = \{^j * fe (l,2,2)(y), 



V n(x^...i%h n }(v) = 

10 V 5 



1 /3\ m 
2 9(xm,s')W{y) = -(-) * ro (3,l,l)(y). 



Substituting these equations into the most recent equation for G(x m ,y) completes 
the proof. □ 



Lemma 8.4. The Green's function evaluated at z m is 

(8.6) G(z m ,y) = -(-) ^f t (l,2,%) + - -)^3^(0,-l,l)(i/). 



10 V5/ ^ 10 V 5" , 

v 7 fe=i v / fe=i 

Proof. We use a similar process to find the formula for G(z m ,y). Note the following 
observations: 

(1) If fc > to, then tp% (z m ) = 0. 

(2) If fc = to, then ip Zm (z m ) = 1. If fc = to and s 7^ z m , then ip™(z m ) = 0. 

(3) If fc < to with s ^ yk and s 7^ Zfc, then ^J(z m ) = 0. 

Using these facts, we have 

m— 1 
G(Zm,y)=E X! [9(Vk,s')ll>y k {Zm)+g{Zk,s')ll> Zk {z m )]lps>(y) 

fc=i s 'ev fc \v fc _i 
+ 2 <Kw'W(v)- 

s'<£V m \V m -! 
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Using the harmonic algorithm, for k < m, we have 

-. /o\ m ~ k 1 / 1 \ m ~k -. /o\ rn ~ ^ i / i \ 'tn — k 

^^=2U) 2 U) ^ ^ (Zm)= 2UJ + 2U 

Since <7(s, s') = if s and s' are in different cells of level k — 1, we deduce that 

fe 

10 V 5, 



J2 9(yk,s')iP k Ay) = ^(l) * fc (l,3,l)(i/), 

4\V fc -i V ' 

]T fl(^y)^(l/) = ^f|) **(1,1, 3)(y), 



10 V 5, 

1 /"}\ m 
£ ^wO^M^U) * m (l,l,3)(|/). 

s'ev m \v m _i ^ ' 

Making these substitutions into the previous equation for G(z m ,y) completes the 
proof. □ 

Lemma 8.5. If u = ouVq and Am exists on SG, then 

(8-7) 

cU(x m ) - § ( -^ J S 3 fc / * fc (0, -1, 1)A« dy-\ [ * TO (1, -1, 1)A« d V -¥J m , 

\ J / fc=1 JSG 2 JSG 

where tp m = /„ ip Xm Au dy. 

Proof. Let w be the harmonic function on Z m determined by the boundary values 
v(x m ) = 1 and v{z m —i) = ^(zm) = 0. Note that v = ip Xm on Zm- Since Z m is a cell 
of level to and v is harmonic, using (j4.ll) with the proper normalization constant, we 
have <- d n v{x m ) = 2(5/3) m and \ d n v{z m _ 1 ) = \ d n v(z m ) = -(5/3)™. These 
equations, together with the symmetric Gauss-Green formula (I4.3[) applied to the 
functions u and v, yield 

•s- d n u(x m ) = / V^Audy+l-j [2«(x m ) - u(z m ) - u(z m _i)] . 

Using the Green's formula, the formulas for G(x m ,y) and G(z m ,y) given by (j8.5|l 
and (J8.6I) respectively, and the normal derivative matching condition at x m yields 
the desired formula. □ 

8.2. Lemmas for Sequences. 

Lemma 8.6. Given a sequence {a m }, ||5 m (5a m +i — 3a m )|J£oc < 00 if and only if 
a m = A(3/5) m + a' m with ||5 m a^J|^°o < 00. Furthermore, 

||5"XJ|,oo <||5 m (5a m+1 -3a m )||^. 

Note that the equation for a m and the bound for a' m implies A = lim m ^ 00 (5/3) m a m . 

Proof. Clearly the second statement implies the first statement. Conversely, making 
the substitution d m — (5/3) m a m , we find that 

3||3 m (d m+ i - d m )\\£oc = ||5 m (5a m+ i - 3a m )||£oc < 00. 
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This inequality implies that {d m } is a Cauchy sequence and by completeness of the 
reals, d m -> D for some 25. Then a m = (3/5) m 25 + (3/5) m (d m - 25). Writing d m 
as a telescoping series 



d m = D + y^ (dk - rffe+i) 
and using the inequality ||3 m (rf m+ i — <2 m )||^°° < oo, we obtain 

oo 

Mm - D| < 2 |d fc - dfc + i| < ^||5 m (5a m+ i - 3a m )||^ 

Then defining a^ = (3/5) m (d m — 25), we see that 

||5"XJ|/oo = \\3 m (d m - DJII/co < ||5 m (5a m+1 - 3a m )||*> 



D 



Lemma 8.7. Given a sequence {rj m }, ||3 m (3ry TO +i — ?7m)||£cc < oo if and only if 

P m Vm\\ Llp <oo. In fact, 

p m (3v m+ l-Vrn)\\^=Wi m Vrn\\ Llp . 

Proof If ||3 m (37? m +i - i] m )\\ eoo < oo, then 

P m7 7m|lLip = SUp3 m |377 m+ i - T) m \ = ||3 m (3?7 m+ i - f]m)\\eoo < OO. 
m 

Conversely, if ||3 m 77 m || Li < oo, then 

-. . . _ q™^ I <-- iiq™^ I 

Lip 



3 m |3r7 m+1 - 3t7| = |3 m+1 7? ro+ i - 3 m Vm \ < ||3 m r; m || L . < 00. 



D 

8.3. Lemmas for Series. 

Lemma 8.8. Fix a constant r < 1 and a sequence {a m }. Then \\r m ' 2 a m \\p < 00 
if and only if ||r m ' 2 (a m +i — a TO )||/2 < 00. More specifically, 

\\r m/2 a m \\ £ 2 < d| ai | 2 + C 2 \\r m / 2 (a m+1 - a m )|| £2 . 

Proof. The first statement obviously implies the second statement. Conversely, 
writing a m as a telescoping series 

m — l m— 1 

a m = a\ + 2_^ ( a k+l — dk) = Ol + 2^ (flm-fc+1 — Om-fe)) 
fc=l fe=l 

we see that 



- :->„ = r ™/2„, a. ^ ( flm _ fc+1 - a m _ fe )r( m - fc )/ 2 r fe / 2 . 



r ' a m = r ' 01 + 

fe=i 



30 
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Using Minkowski's inequality, we have 

m— 1 

E( a ™-*+i- a ™-^ (m ~ fe)/2 r fe/2 
fe=i 



e 



- X! rfe/2 ( a m-k+i - a m -k)r 



(m-fe)/2 



Xfc<r 



fc=l 



<]Tr fc/2 ||(a„ l+ i -a m K»/ 2 1|, 2 . 



fc=i 



Using Minkowski's inequality again and the above inequality we find that 



\\r m/2 a m \\ l ,<\\r m ' 2 a 1 \\ l 2 + 



which completes the proof. 



£(a m _ fe+1 -a m _ fc )r( m - fc )/ 2 r fc / 2 



fc=i 



a 



Lemma 8.9. Fza; a constant r > 1 and a sequence {a m }. Then a m = ^4 + a^ wit/i 

||r m ' 2 a^||£2 < oo if and only if \\r m ' 2 {a m +i — a m )\\f2 < oo. In fact, 

ii m/2 / II ^- y^r||_"^/2/ \|| 

Proof. Clearly, the first statement implies the second statement. To prove the 
converse, we first show that {a m } is Cauchy. For m > n, we have 

m— 1 

a m - a„ = ^ (a/b+i - a fe ) r k/2 r - k/2 



and applying Cauchy-Schwarz yields 



1/2 



|«m ~ On | < I y^ (afe+1 - afe) 2 ?'' [^2 ~k 



1/2 



<C\ — . 



ft follows that {a m } is Cauchy and by completeness of the reals, a m — > A for some 
A. Since 

oo oo 

O-m — A = 2_^ { a k — Ofe+l) = ^2^(a>n+fc ^ dm+k+l)) 



fe=0 



we see that 



r m ' 2 (a m - A) = J2 r(m+k)/2r ~ k/2 ( a ™+* - "rn+k+l). 



fe=0 



Using this equation and Minkowski's inequality, we have 



^ 2 (a m - A)y <J2r- k/2 \\(a k+m - a k+m+1 )r^ +m ^ 2 \ 



c* 



fc=0 



< J2 r k/2 \\( a m-a m +i[ 



r m/2 \U>, 



fc=0 



which completes the proof. 



□ 
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Lemma 8.10. Given a sequence {a m }, ||(25/3) m ' 2 (5a m +2 — 8a m +i + 3a m )\\£2 < oo 
if and only if a m = A x + A 2 (3/5) m + a' m with ||(25/3) m / 2 <J|^ < oo. More 
specifically, 

\\(25/3) m / 2 a , m \\ e 2<C\\(25/3r(5a m+2 -8a m+1 + 3a m )\\ P . 

Note that the equation for a m and the bound for a' m imply that A\ = linim-^ a m 
and A 2 = lim m ^ DO (5/3)" l (a m - A{). 

Proof. Clearly the second statement implies the first statement. To prove the con- 
verse, we apply Lemma I8U1 twice . Making the substitution 3 m d m = 5 m (a m+ i — a m ) 
yields 

oo / ok\ m °° 

5^ ( T ) ( 5flm + 2 ~ 8a ^+l + 3a m) 2 = 9 ^^ 3 m (c?m+l - dm) 2 < OO. 
TTi—1 ^ ' m— 1 

The hypotheses of the lemma are satisfied for {dm}, so we have d m = D 4- d f m with 

oo oo 

^3 m |< n | 2 <C^3 m (d m+1 -d m ) 2 . 

771—1 m— 1 

In order to apply the lemma again, define e m = a m + (5/2)(3/5) m D so that 

oo oo / (-. _ \ m 

]T3™KJ 2 =E(f) (e m+1 -e m f<^. 

m—l m— 1 ^ ^ 

Using the lemma again, except on the sequence {e m }, we have e m = E + e' m with 
the estimate 

Wf) KJ 2 <cWf) {e m+ i-e m f. 

m— 1 ^ m—l ^ ^ 

Finally, using the definition of e m , we find that a m = E — (5/2)(3/5) m D + e' m . 
Combining the above equations and inequalities, we obtain 

oo , s m oo , _. m 

E(f) l e ™| 2 ^ C E(f) (5a m+2 -8 flm+1 + 3a m ) 2 . 

D 

Lemma 8.11. Given a sequence {r] m }, \\3 m ^ 2 (3rj m +2 — 16r/„,+i + 577 m )||£2 < oo if 
and only if i] m = 5 m A + 77^ wii/i ||3 m ' 2 ??m||f2 < 00. Furthermore, 

P m/2 VL\\1 < Ci(7?2 - 5T7i) 2 + C 2 ||3 m/2 (377 m+2 - 1677^+1 + 5ry m )|| 2 2 . 

Proof. The second statement obviously implies the first statement. To prove the 
converse, we use both Lemma 18.81 and Lemma 18.91 Define e m = 3 m (r) m +i — 5?y m ) 
so that 

00 00 _. 

E 3m ( 37 ?™+2 - 16?7 m+ i + 5?7 m ) 2 = J^ ^r( e ™+i - e m ) 2 < 00. 

m—l m—l 

Applying Lemma 18.81 to the sequence {e m } gives us 

OO -. oo 

E ^rM 2 ^ C il e i| 2 + C 2 E 3m ( 3 ^»+2 - 1677 ro+ i + 5r/ m ) 2 < oo. 

771—1 771—1 
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Making the substitution 5 m d m — rj m , we see that 

oo _. oo oo 

E ^rl e ™! 2 = E 3 " l (^+i - 5 V,nf = 25 Y, 75 m (d m+1 - d rn f < oo. 

m— 1 ?7T-=1 m— 1 

Applying Lemma [5TU1 to the sequence {d m }, we find that d m = £> + d' n with 



^ 75 m |<J 2 < C Y, 75 m (d m+ i - d m ) 2 . 

777 = 1 

It follows from the definition of d m that rj m — 5 m D + 5 m d / m . Defining 7f m = 5 m d^ 
and combining the above equations and inequalities, we obtain 

oo oo 

Y, yn WJ 2 < Ci(v2 - 5 m f + C 2 Y 3 m (3r/ m+2 - 16r ?m+1 + 5r, m ) 2 . 

771—1 771 — 1 

□ 
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